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Abstract 


Some models have extremely low-mass pseudo— Goldstone bosons that can 
lead to vacuum phase transitions at late times, after the decoupling of 
the microwave background. This can generate structure formation at 
redshifts z ft 10 on mass scales as large as M ~ 10 18 M©. Such low 
energy transitions can lead to large but phenomenologically acceptable 
density in homogeneities in “soft topological defects” (e.g., domain walls) 
with minimal variations in the microwave anisotropy, as small as ST /T £ 

10 -6 . This mechanism is independent of the existence of hot, cold, or 
baryonic dark matter. It is a novel alternative to both cosmic string 
and to inflationary quantum fluctuations as the origin of structure in the 
Universe. 
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The problem of generating cosmological structure (galaxies, clusters, voids, pecu- 
liar velocities, etc.) in a universe that appears very homogeneous and isotropic on 
the largest scales is, perhaps, the major problem in physical cosmology today. Pre- 
viously, attempts at solutions have involved generating density fluctuations (either 
quantum mechanically, the usual Gaussian inflationary fluctuations; or topologically, 
such as cosmic strings) at a very early cosmological epoch (e.g., Grand Unified [GUT] 
epoch when kT ~ 10 14 GeV) which survive to serve as seeds at the galaxy formation 
epoch at kT 10 -2 eV. In some scenarios thes e seeds gravitationally accrete large 
quantities of non-baryonic dark matter, whereas in others they explode and push the 
baryons about. Rather than go into a detailed commentary on each of these models, 
let us merely note that the following combination of observations has been difficult 
(but maybe not impossible) for any existing model to satisfy: 

1. Microwave background anisotropy 1 ST /T 2 x 10 -5 . 

2. Quasars and some galaxies exist 2 by redskifts 2^3. 

3. Large scale peculiar velocities 3 depart from the Hubble flow up to Av ~ 600 km/ sec 
on scales of R ^ 30 Mpc. 

4. Structure, clustering, foam, voids, etc. exist 4 on scales of ^ 20 Mpc. 

5. Clusters of galaxies appear more strongly correlated 6 than galaxies. 

People favoring some theoretical models may dismiss some items so as not to have to 
discard their model prematurely. For example, proponents of cold dark matter and 
Gaussian fluctuations don’t like #3 and #5, and most Gaussian fluctuation models 
with linear growth must argue that $2 is a statistical tail and not representative of 
the bulk of galaxy formation. Moreover, #1 forces a fine-tuning of the GUT phase 
transition. 

The purpose of this note is to propose a completely different alternative, where 
the fluctuations are generated after decoupling, thus inducing minimal fluctuations in 
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6T/T. Furthermore, the fluctuations will be associated with “soft” topological struc- 
tures, having fixed internal energy densities typically of order m^ u(rtrio . At some 
redshift z<Ucoupiing > zi > 1 we will have 6pf p ~ 1, so density fluctuations immedi- 
ately grow nonlinearly subsequent to Z\. This avoids the long linear growth period in 
previous models that makes it relatively difficult to have galaxies form from small ini- 
tial amplitude by redshifts z 1. Large variations in density can also produce large 
peculiar velocities on scales up to the horizon at the time of the transition, ~ 100 Mpc 
today (in present distances, the horizon at redshift z 1 is 3000z -1 / 2 Mpc). Observed 
structures as large as 1O 18 M 0 have been claimed by Tully 8 and cannot be easily ac- 
commodated by models where fluctuation scales must be limited by the horizon scale 
at decoupling. Another possible feature of our mechanism is that the fluctuations 
associated with the critical point phenomena of a phase transition have a scale-free 
or fractal character out to the horizon at the epoch they form. Such a fractal seems 
to provide a natural way of understanding the cluster-cluster correlations 7 as well as 
producing patterns that resemble the observed large scale structure. Much of what 
we say here is generic to any late-time phase transition (including, e.g., one in the 
hidden E% sector of E 6 x E% superstring theories that could create the observed struc- 
ture and be otherwise unobservable), though in the Appendix we present a specific 
model, intended primarily as an “existence proof” of the viability of such schemes. 

Our initial thought as to why one might expect such a late-time phase transi- 
tion came from considerations of the possible cosmological implications of low-mass 
pseudo-Goldstone bosons which can arise naturally in a variety of GUT settings. 
Pseudo-Goldstones, such as massless familons 8 , arise when the pattern of masses of 
the observed fermions is associated with a spontaneously broken, continuous (but 
ungauged), exact symmetry. With further small explicit breakings of these symme- 
tries, familons acquire minuscule masses, e.g., in the “schizon” models of Hill and 
Ross® these are typically of order ~ mj//*, where /* ~ 10 15 GeV to 10 18 GeV 
is a generic GUT scale, and m/ the mass of the associated family of fermions. The 
resulting masses of these pseudo-Goldstone bosons can have interesting astrophysi- 
cal implications if the associated fermion f ami ly is taken to be the neutrinos. What 
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are the neutrino masses? An examination of MSW 10 mixing solution to the solar 
neutrino problem suggests neutrinos have a small mass, m v ~ 10~ a eV With these 
typical values for neutrino masses, one thus estimates the Compton wavelength of the 
neutrino— schizon to be r = ~ 100 pc to lMpc (for 10 1 to 10 eV neutrino 

masses), with /* ~ 10 1B GeV. This is a remarkable result for a cosmological distance 
scale, coming from the ratio of two high-energy particle physics mass scales and is 
suggestive of a possible mechanism for formation of large-scale structure. 

What dynamics can lead to a late time phase transition? In the Appendix we 
write down the low-energy effective Lagrangian for a particular scheme and find that 
the potential for the neutrino-schizon, (f>, is given by: 

V{<j>) = -c(T)m*cos(2<f>/f) (1) 

which implies m\ = 2c(0)77i*//£. Here c(T) is a temperature dependent coefficient 

of the form ~ log(T//x)/47r 2 , valid for T Z m v and c(0) log (fi/m^/ATT 2 , where 

pi is a renormalization scale as described in the Appendix. We shall assume that 
pi, Z mv and that the initial state of the <f >~ field is given by a quantum-mechanical 
vacuum (a gaussian wave-functional in <f>) localized around a field configuration in 
a minimum of e.g., (<fi) = 0. We emphasize that this need not require that 

the pseudo-Goldstones are in thermal equilibrium (they generally are not); the other 
matter fields contribute to c(T) through Feynman loops with thermal expectation 
values so long as they have a thermal density matrix (they themselves need not be in 
equilibrium; the thermal density matrix applies after decoupling). 

In general as the universe cools, the pseudo— Goldstone bosons can undergo a phase 
transition if, for example, at some temperature T c the coefficient c(T) can become 
zero and, for subsequent lower temperatures, change sign. Considering the particular 
model described above this phase transition occurs at T e = /i. We define this to cor- 
respond to a redshift z 0 , and we shall further assume that c(T) evolves monotonically, 
as is the case above. The pseudo— Goldstone field now acquires a new VEV, with a 
discrete degeneracy that leads to topological defects, here typically domain walls. It 
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is not unnatural to choose fi ~ A m u , with A a dimensionless parameter that can range 
as small as C?(l); we shall typically consider A ~ 10. A prediction for fi depends upon 
the precise details of the model at the GUT energy scales. The fact that the late 
time phase transition can, in principle, naturally occur is sufficient for our present 
discussion. 


This phase transition most likely proceeds quantum mechanically, in close analogy 
to inflationary phase transitions. The wave-function, to find a given Fourier 

component of the pseudo-Goldstone field, <f>k, satisfies the Schroedinger equation: 


1 

2 


S 2 

6*1 


+ {k 2 + HO M 


*(<M = 


( 2 ) 


(where we assume localization around 0* = 0 in writing the F"(0) term; this term 
is time dependent through c(T) and T(t)). As the potential V(<j>) changes sign the 
wave-function uniformly spreads in <f>k provided the momentum k is not too large. 
Using the potential of eq.(l) we have spreading provided, 


k 2 - 4|c(0)| mj/ij < 0 


( 3 ) 


Of course, some components corresponding to small k will begin to spread as T ~ 
T c , but ultimately as T — ► 0 all will spread with k satisfying eq.(3); this 

relation establishes a critical kf = 4c(0)mJ//^ above which no spreading occurs. 
Spreading of the wave-function continues for a time t ~ /^/m* until it becomes 
equally likely to find a new VEV, either v\ = tt/^ or v a = — irf^, Thus, some regions 
of space will have the wave-function collapse around t>i, while others around vj. There 
will necessarily be domain walls between these regions. The domain walls have the 
transverse structure of the kink— soliton of the sine— Gordon equation: 


-T- = 4tan -1 


exp (m,p(x - /3t)/J[l -/?*)) - ir 


( 4 ) 


with v/c = fi. The thickness, r, is thus ~ /^/c(0)m* which is also of order 
the initial characteristic spacing at the time of the phase transition, h' 1 and is the 
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characteristic time scale for the duration of the phase transition. The thickness of 
these soft domain walls can range from 100 pc to 1 Mpc for 1 to 10 eV neutrino 
masses, with / ~ 10 15 GeV. Since the domain wall thickness is of order A" 1 at the 
initial redshift of the phase transition, z 0 , the field configurations after this transition 
are very ill-defined and the domain walls are hard to distinguish and are distributed 
randomly, often lying on top of one another. The quantum state remains highly 
excited. It will relax, however, by Hubble red-shifting and the wall interspacing 
will typically become fc c _1 (l + *o)/( 1 + -*) at subsequent redshift, z < z 0 , while the 
walls themselves will relax to the well-defined kink-soliton configurations of eq.(4). 
In addition, there may be slow recombination of structures and evolutionary effects, 
analogous to the evolution of cosmic string networks. 

At the epoch of the phase transition, a redshift z 0 , the average cosmological energy 
density in the <£ field configurations is of order <0)mj. This energy density will now 
redshift, the leading contribution behaving like a domain wall which we write as: 


Pwallt ~ cl 




l + K 


( 5 ) 


where the parameter k > 0 is introduced to model the various evolutionary effects 
(parenting processes, domain wall annihilation and decay, etc.; there are, of course, 
subleading contributions to eq.(4), e.g., the <fi field will oscillate about the true vacuum 
value between the domain walls, in analogy to the axion field, but this energy redshifts 
as (1 + z) 3 and may be neglected). On the other hand, matter redshifts as 


P matter — Pmatter{ z l) ^ ^ J ' ' 

In the center of a domain wall we would have a fixed, non— redshifting energy 
density of order c(0)m*. Let us assume that at the special redshift, z i, the matter 
energy density becomes equal to c(0)m* (this is determined once the value of the 
neutrino mass is specified, thus for m u ~ 10“ 2 we have z\ ~ 10). Indeed, Z\ can in 
principle equal z 0 , but we expect in general that Z\ < Zo- At the redshift z\ the domain 
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walls will appear as 0( 1) density contrasts, Sp/ p ~ 1 and large-scale structure will 
begin to form rapidly. Comparing eq.(5) and eq.(6) at z = 0 gives a relationship for 
the closure fraction of domain walls to matter: 

ttwolli _ (1 + 3j) 3 

^matter (id" 2 o) 1+ * 

This implies that we have a closure constraint in the universe of the domain wall 
energy density, £ 1, with fi motter « 0.03 to 0.10, implying that the redshifts z 0 

and z\ satisfy: 


30 to 10 £ z\jz^ +K) 


( 8 ) 


For example, in the extreme case that we close the universe with domain wall energy 
by saturating inequality (8), and further assume k js 0 we might expect zq ~ 10 2 , 
and then z ^ ~ 10 is the redshift at which Spjp ~ 1, which requires that m v ~ 
(pmatterizi)) 1 ' 4 ~ lO" 2 eV. 

The important point here is that, because the fluctuations are immediately non- 
linear, at 2 ~ zi, they will immediately grow into larger structures. The scenario 
is viable regardless of whether or not dark matter is present, or if the dark matter 
is of the hot or cold variety. Moreover, by having a phase transition occuring after 
microwave background decoupling at kT ~ leV, current STjT limits are not so 
restrictive. 

However, ST /T is roughly related to /* as follows. Consider the shift in energy 
of light falling into the center of a static domain wall with central energy density 
p ~ c(0)m* of characteristic thickness 7* rsj /^,/c(0)m 2 . We thus have: 


ST 

T 


Gpr 2 ~ 


fl 

Mpianck 


( 9 ) 


This is actually a crude upper bound for the net ST/T generated in this scenario, 
since there will generally be a cancellation of blue— shift against red— shift in the limit 
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of a completely adiabatic transition. For us, the value of is, therefore, probed by 
measurements of 8T/T while 8p/ p is independent of f^. Thus, in the present scenario 
8p/p can be large at z i, while 8p/p in other scenarios, such as cosmic strings, is 
constrained by the upper limit on 8T/T to be small. 

Indeed, conventional phase transitions occuring prior to decoupling are plagued 
by the requirement that all fluctuations be small. 11 Hence, associated defects like 
domain walls were required to be outside of the horizon or else to be “inflated” 
away. The inverse of this problem is that other topological structures like cosmic 
strings which also form at GUT scales would require fine tuning to avoid also being 
inflated away (only strings forming after inflation survive — but why should a GUT 
string form after inflation?). A late-time phase transition producing soft topological 
defects, however, has the advantage of producing structure without a conflict with 
inflation. The phase transition for any such low mass particles will be analogous to 
the axion phase transition that has been well studied theoretically. Wasserman 12 , in 
a preliminary treatment of late time phase transitions without an underlying particle- 
physics motivation focused on the small perturbative fluctuations. We focus instead 
on the large topological defect effects; a real late time phase transition would have 
both, however the topological defects have the advantage that a fixed non-redshifting 
energy density is available in the cores of the defects which allow 8pj p ~ 1. The energy 
scale of the topological defect is relatively low; hence the structures are “soft,” e.g., 
their core sizes are very large and they have internal energy densities of order m * in 
the schizon model (compared to GUT strings or monopoles having microscopic cores 
with energy scales of ~ (10 16 ) GeV). Nonetheless, the density contrasts can be high, 
after some redshift z x , with domain walls or membranes having variations relative to 
the surrounding true vacuum of order unity or greater. 

How the generic soft topological defects actually evolve, thus providing an esti- 
mate of k in the above expressions, is a difficult question and no doubt varies with 
model specifics. As in the case of cosmic strings, the evolution of domain wall or 
“cosmic membrane” networks is non-trivial (in the string case, much debate still 
continues on this problem). The soft cosmic membranes can both accrete or repel 
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matter. 14 To nearby matter, because of large surface tension, a domain wall is grav- 
itationally repulsive 14 ; this will mimic anti-biased, hot dark matter. 17 Closed bubbles 
have positive mass as seen from outside at distances greater than the radii of the bub- 
bles, and will serve as accretion centers just as loops of cosmic string. If accreting, 
then the subsequent evolution of the gravitating baryons and dark matter is that of 
fragmenting sheets, but with collapse times that are much more rapid than standard 
linear growth. If repelling, then the matter between the walls will be compressed and 
collapse, somewhat in analogy to the scenario of Ostriker et al. 1B In fact, repelling 
late-time cosmic membranes may enable an explosive galaxy formation mechanism 
to work without having to invoke primordial magnetic fields and superconducting 
strings, or other contrivances. The initial pattern of membranes will be determined 
by the zeroes of a random field, and thus will resemble the pattern of caustics seen in 
hot dark matter models 16 . Finally, scale-free critical fluctuations within one domain 
may provide yet another possible source of large-scale structure. 

As mentioned above, the microwave background anisotropies expected in this class 
of model are remarkably small. To first order they are zero, since no fluctuations need 
exist prior to electromagnetic decoupling at the time of recombination; in fact, the 
surface of last scattering for the microwave photons could be perfectly smooth. (This 
is even less than with cosmic strings, where string induced gravitational fluctuations 
exist at last scattering.) The only ST jT effects that are induced are due to the second 
order differential blueshift/redshift resulting from the propagation of the microwave 
photon through evolving transparent density fluctuations. As shown in eq.(9) the 
fluctuations induced by the domain walls are limited by roughly fl/Mj> lanck . Present 
observational limits on ST jT thus constrain ^ 5 x 10 16 GeV which clearly allows 

to be of GUT scale. In fact, note that for 10 16 GeV the maximum implied ST jT 

due to the domain walls themselves is 10~°. Once we consider such small ST jT 
we must be careful not to ignore the effects due to the eventual propagation of the 
microwave photons through the transparent potential wells of the created structures. 
While we need not consider the Sachs- Wolf effect on the surface of last scattering, 
Rees and Sciama 18 note that the present existence of large structures will produce 
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a second-order effect in 8T/T for fl = 1 universes. However, the magnitude of this 
effect we estimate to be only 8T/T ~ 6 x 10 -8 Mi 8 hgiZj oo where M 18 is the mass 
of the structure in units of 1O 18 M 0 , Rioo is the size of the structure in units of 
100 Mpc, and h 0 is the usual Hubble constant in units of 100 km/sec/Mpc. Thus, 
even the largest structures claimed by the most ambitious observers 6 would not yield 
observable effects. Note that even if fl ^ 1 or induced Sachs-Wolf effects enter^late- 
time transitions would still only yield effects of order 8T/T ~ SM/Mo where 8M is 
the mass of the fluctuation and M 0 « lO 33 ^ 1 Af 0 is the horizon mass today. Thus, 
the largest scales presently observed where 6M ~ lO 17 Af 0 only yield 8T/T ~ 10 -6 , 
and then only on angular scales Ad ~ (a few degrees), or larger. All smaller scales 
would yield even smaller effects. It is clear that a late time phase transition produces 
the smallest possible 8T/T of any proposed scenario. If measurements of 8T /T are 
ever reported at levels of a few xlO -6 , they would directly constrain and possibly 
measure via eq.(9). 

Observationally there is also a limit on the fraction of critical density in the <f> 
field today, due to the induced large scale velocities. From the present data 3 on 
R ~ 40 Mpc we know that fl^(l + z$)8pjp £ 0.2, where z $ is the redshift of the phase 
transition and 8pj p is the density variation in the $ field. Thus, for 8p/ p ~ 1 we 
have < 0.2/(l + z#). This constraint sets bounds on the evolution of the <f> field 
structures, including k from eq.(5). 

The recent observations 19 of an excess at sub-millimeter wavelengths in the mi- 
crowave background may, if real, also be explained with the help of a late-time phase 
transition. In particular, this non-linear growth model may be the only way to have 
significant star formation at z ^ 30. Hogan, et.al. 30 argue that such star formation 
could create the necessary ionization. It should also be noted that energy released 
by the phase transition itself or by decay or annihilation of topological defects might 
provide an alternate source for ionization. 

Laboratory tests for the model vary with the specific details. If the late-time 
transition is associated with neutrino-schizons then there must exist small neutrino 
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masses near the MSW range. The GALLEX and Baksan gallium experiments coupled 
with the 3r Cl experiment and the proposed DjO experiment will eventually confirm 
or deny the MSW explanation for the solar neutrino puzzle. If MSW, then are there 
neutrino associated schizons? Our mechanism requires a generic pseudo-Goldstone 
boson which will be hud to detect directly, but its brethren associated with charged 
leptons or quarks produce potentially observable new phenomena, e.g., new composi- 
tion dependent pseudo-Gravitational forces, as detailed in ref. (9). The observation 
of such effects and non-zero neutrino masses would be compelling circumstantial ev- 
idence for possible cosmological effects proposed here. As a generic mechanism, it 
could even occur in the hidden sector of JE7g x E& superstring theories and be impos- 
sible to observe except by its role on galaxy formation. 

Obviously much work remains to be done to examine the details of this class of 
models. In particular, the astrophysics of the detailed large scale structure that is 
generated by such late time fluctuations is only sketched here; and full hydrodynamic 
calculations will have to be carried out. Furthermore, detailed particle physics models 
will have to be developed to see if all the preferred properties really exist in a fully 
consistent model. Eventually we would hope to make detailed quantitative predictions 
about the model vis-a-vis large scale structure. However, the present large scale 
structure observations are still quite qualitative. Quantitative statistical measures 
have yet to definitively describe the apparent structure in & reproducible manner. 
Anecdotally, voids, filaments, sheets, bubbles or sponges appear, depending on the 
analyses used and on the rapporteur. Conceivably cosmic membranes could make 
any or all of these structures depending on how they evolve. Hopefully, specific 
quantitative predictions will be made before the observational data converge. Our 
purpose here is to alert readers to the fact that an alternative to the standard galaxy 
formation scenarios may exist. The physics it relieB upon is not any more exotic than 
the GUT physics that the standard scenarios utilize. At low energy scales the model 
might even be testable in the laboratory. In any case, it may be the only model that 
can survive limits on ST jT & 10 -8 . 
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Appendix 

Here we present a simple neutrino-schizon model which can produce a phase 
transition at late times in a natural way. We do not discuss the fundamental Higgs 
structure at the GUT level which allows such a model to be natural; this can be done 
along the lines of the discussion in ref.(9). 

We assume two species of neutrino, i/\ and and a pseudo— Goldstone boson 
£ with decay constant f $ - <f> appears as a phase factor in the low energy effective 
Lagrangian. The mass terms for the neutrinos are assumed to take the form: 

i'lL^iR (w* + ee**^*) + vilVir ( m — + h.c. (A.l) 

The special choice of eq.(A.l) actually corresponds to a discrete symmetry under 
which 1/1 <-► v 2 and <f>/ -*• <£//* + 7r. In the limit m = 0 the theory has a continuous 
chiral symmetry realized nonlinearly with <f> in analogy with the Peccei-Quinn sym- 
metry and the axion (actually, we may associate <f> with the m terms by a redefinition 
of the neutrino fields, so the chiral symmetry may be regarded as present in the limit 
e = 0 as well). The presence of both nonzero m and e implies the explicit breaking 
of the chiral symmetry; the <f> field will acquire a mass. Note that we can rewrite 
eq.(A.l), after a ^-dependent phase redefinition of the v as: 

V\V\ (m 3 + e 3 + 2 me cos <f >/ /*) 1/3 + P 2 u 2 (ru 3 + e a - 2me cos <j>j /*) ^ ( A.2) 

The discrete symmetry protects the mass of the <f> field from being quadratically 
divergent when loops are considered. However, we do obtain an induced logarithmi- 
cally divergent mass term for <f> of the form: 

— m 2 e 2 ln(A/m) cos(2<£//^)/47r 3 (A. 3) 

Thus, we must introduce a renormalization coimterterm for the induced <f> mass. In 
principle this can be arbitrarily large, but it is not an unnatural fine-tuning to choose 
a counterterm of order the result of eq.(A.3); viewed another way, eq.(A.3) implies 
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a renormalization group equation for the <f> mass such that, if we define the mass to 
be zero at some large energy scale A, it will never grow large at low energies and will 
be given more or less by eq.(A.3). We thus have for the renormalized 4> mass the 
expression: 

— m 3 e 3 ln(/i/m) cos(2^//^)/4ir 3 (A. 4) 

where n summarizes the choice of renormalization condition. 

We now wish to examine the system at finite temperature (Note: this does not 
imply that the neutrinos are in thermal equilibrium; it merely implies that they exist 
at finite density in a thermal density matrix, which can be a relic of an earlier epoch 
in which they were in equilibrium). The neutrino bilinears at high temperature have 
the behavior: 


( vv) T —* mJT 2 / 4tt 3 + m 3 ln(T/m 1/ )/87r 3 T > m„ (A. 5) 

Here m„ is the full physical mass of the neutrino as read off in eq.(A.2): 

m„ i = (m 2 + e 3 + 2mecos f 4 ) * ; m v3 = (m 2 + e 3 — 2mecos ^ (A. 6) 

Thus, substituting the finite temperature expectations of eq.(A.5) into eq.(A.2) and 
adding the zero-temperature <p mass of eq.(A.4), gives the temperature dependent <f> 
mass: 


m 3 e 3 ln(T//i)cos(2^//0)/47r 3 (A. 7) 

where the T 3 terms have cancelled (these are analogues of the quadratic divergence 
at zero temperature and cancel owing to the discrete symmetry; note that the overall 
sign of eq.(A.7) is irrelevant as we can always shift — ♦ <f >/ + tt/ 2). 

We see from this result that the temperature of the phase transition is controlled 
by the renormalization mass, fi. For T > fi > m v the potential has the positive sign, 


-13- 


FERMILAB-Pub-88/120-A 


and for m v < T < p, the negative sign. If we choose e ~ vn v ~ 10 2 eV, /* ~ 10 1B 
GeV and p ~ 10m„, we have the scenario outlined in the text of this paper. 

The unsatisfying aspect of this is that we do not have a prediction for the quantity 
p. This would have to come from a detailed understanding of the full GUT theory, 
which we do not know. There is nothing in principle wrong with a value of p ~ m v , 
though one might say that we are making a special choice of the log-interval of p ; 
however, this is not too special a choice since p between 10 2 to 1 eV is acceptable, 
and this is therefore only a particular selection of 1 in 13 log— intervals between 10 
eV and 10 24 eV, the GUT scale. Nevertheless, this model (taken together with a 
demonstration of the naturalness in the full GUI theory as in ref.(9)) illustrates that 
a late-time phase transition is not at all unreasonable and may even be dictated in 
some models. The above model is simply a toy; there are no doubt large classes of 
models admitting this phenomenon. 
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